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We present an experimental and theoretical study of the magnetic field dependence of the critical 
current of Josephson junction ladders. At variance with the well-known case of a one-dimensional 
(ID) parallel array of Josephson junctions the magnetic field patterns display a single minimum 
even for very low values of the self- inductance parameter /3l. Experiments performed changing 
both the geometrical value of the inductance and the critical current of the junctions show a good 
agreement with numerical simulations. We argue that the observed magnetic field patterns are due 
to a peculiar mapping between the isotropic Josephson ladder and the ID parallel array with the 
self-inductance parameter /3£ ff = /3l + 2. 



74.50.+r, 85.25.Dq 



I. INTRODUCTION 

A Josephson junction ladder is an array of coupled 
superconducting loops containing small Josephson junc- 
tions as shown schematically in Fig. |l](a). In the past 
years, ladders of Josephson junctions have attracted con- 
siderable interest for a number of reasons. On one hand, 
more complex systems such as two-dimensional arrays of 
Josephson junctions can b£_viewed as elementary lad- 
ders coupled to each otheitra. On the other hand, the 
ladders are an ideal ground for experimental and-theoret- 
ical investigations of discrete nonlinear entitiestHj, such 
as breathersEliiy and vortex propagation^. In spite var- 
ious groups have numerically and theoretically studied 
the static properties of Josephson ladderstlfEj, so far no 
systematic comparison with experimental data has been 
carried out. Grimaldi et al.E3 have performed numeri- 
cal simulations on these systems. Their findings are that 
the behavior of ladders is quite different from that of 
one-dimensional (ID) parallel arrays. In contrast to a 
ladder, a ID parallel array contains only Josephson junc- 
tions in the direction of the bias current Jb but no one 
transverse to it. A ID parallel array of Josephson junc- 
tions placed in magnetic field shows a pattern of critical 
current with as many minima as the number of loops 
(for relatively low /3l, as defined below)Ej'L3. In Joseph- 
son ladders with junctions on the horizontal branches not 
only the number of minima in the pattern does not corre- 
spond to the number of loops even for extremely low /3l, 
but also the pattern dependence on the parameter /3l is 
different: the critical current Iq never reduces to zero for 
fully frustrated arrays (i.e. when there is half flux quan- 
tum in each cell). Baharona et al.U have shown that one 
can analytically estimate the depinning current of fluxon 
trapped into the ladder in the limit of zero inductance. 
They have also computed the onset of instability in the 
case of no fluxons, thus retrieving analytically the numer- 
ical result of Ref. fL2 for very low inductance. Moreover, 



the authors of Ref. [IT] have estimated that the critical 
current of a ladder with a fluxon trapped in each second 
cell is higher than the depinning from the empty ground 
state for a moderately high magnetic field. 

The aim of this paper is to present experiments per- 
formed on isotropic Josephson ladders with various values 
of the self-inductance parameter /3l- We call as isotropic 
a ladder consisting of identical junctions. We also make 
an analysis of the model to explain the observed depen- 
dence of the pattern upon The work is organized 
as follows. In Section II we describe a model for the 
Josephson ladders, in Section III we show the experimen- 
tal findings and make the comparison with the numerical 
predictions. Finally, Section IV contains a discussion of 
our results and Section V the conclusion. 




FIG. 1. (a) The electrical scheme of a Josephson junction 
ladder; crosses (x) indicate Josephson junctions, (b) Optical 
image of one of studied samples. 
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II. THE MODEL 

To derive the equations for the ladder, we start from 
the fhrxoid quantization over a cell: 



cell 



2tt 
$0 



indl 



(i) 



mean to find the current point at which the static solution 
becomes unstable. The value of the dissipation used in 
the simulations is fictitious, and it has been chosen equal 
to 1 for computational convenience. The experimentally 
studied arrays are actually underdamped and therefore 
have a much smaller damping coefficient, a ~ 0.005 - 
0.08, that can be controlled by temperature. 



where <p are the phases of all the junctions in the cell, 
$o is the flux quantum, and <I> cxt and $ md are the ap- 
plied and induced flux, respectively. To evaluate <& md 
we retain only self-inductance terms, i.e. we assume that 
<j,md _ ^js^ js bgirrg the screening current circulating 
in the elementary cell and L is the self inductance of 
the cell. For the junctions we assume the RCSJ model 
and suppose that all junctions are identical (isotropic 
ladder). With these ingredients it is possible to derive 
the following set of equations for the gauge invariant 
phase difference across the vertical and horizontal 
(V>i) junctional: 

4>i + a(f>i + sin 4>i = 
= 7T [<t>i-i - Hi + &+1 + 2(Vi - ipi-i)] +7, i — 2, . . . , N 



4>i + aipi + sin ipi = 
= — [<j>i - - 2iPi + 2tt/] , i = 1, . . . , N. 



(2) 



The boundary conditions are: 

0i + a(f>i + sin 4>x = — [</> 2 - 4>i + 2if>i + 2nf] + 7, 

PL 

4>N+\ + a<fi N+ i + sin^jy+i = 
= -r\<t>N- 4>n+i + 2tp N - 2tt/] + 7. (3) 

PL 

Here, 7 = I-q/i^ t is the bias current normalized to the 
single vertical junction critical current, pY = 27rLz™'/<I>o 
is the self-inductance parameter (the self-inductance L 
of a square cell with the side a can be estimatedlij as 
L = 1.25/ioa, where /io is the magnetic permeability), the 
ratio T] = ic° r /*c er between the horizontal and the vertical 
junction critical currents is the anisotropy parameter, a is 
the normalized dissipation parameter, / = $ cxt /$o is the 
normalized external flux often noted as frustration, and 
N is the number of loops. For the static case the parallel 




correspond to the limit 
) we take advantage of the 



arrays considered in Rcf. 
77 — > 00. In deriving Eqs. (| 
fact that, due to the symmetry of the system, the current 
flowing in the top and bottom horizontal junctions of the 
same cell (see Fig. |l|(a)) differs only in direction but not 
in amplitude, and therefore we can write the equation for 
just one of themEl. Finally, we want to stress that in this 
work we are interested only in the transition point from 
the static to the dynamic solutions (the critical current), 
therefore the dynamics is nothing but a computational 



III. EXPERIMENTAL J c VS / PATTERNS 

We present an experimental study of 10-cell Josephson 
junction ladders. Each elementary cell of the ladders con- 
tains 4 identical small Nb/Al-AlO^/Nb Josephson tunnel 
junctionst£l, which have an area of 3 x 3 //m 2 . To get dif- 
ferent values of pY , we used samples with different critical 
current density j c (100 A/cm 2 or 1000 A/cm 2 ) and also 
varied the loop size a (2.8 fim or 9.9 um). A SEM image 
of a typical ladder is shown in Fig. |l](b). 

We have measured the ladder critical current, Iq, ver- 
sus frustration / for 4 isotropic ladders (77 = 1) with dif- 
ferent pL- The values selected are pY = 3, 0.88, 0.25 and 
0.088, a range where the peculiar behavior of the lad- 
ders should be clearly visible. The measurements were 
performed in a cryoperm shield. The magnetic field H 
applied perpendicular to the substrate was provided by 
a coil placed inside the shield. The uniform bias current 
Jb was injected at every node via on-chip resistors. The 
voltage across the first vertical junction was measured 
to define the depinning current of the ladder. Finally, 
the 7c v [S-i/ dependencies were measured using GoldExi 
softwareO. 
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FIG. 2. Experimental (symbols) and numerical (solid lines) 
7 vs / patterns of ladders with different cell sizes. Parameters 
are: N = 11; /3 L = 3.0 (squares), 0.88 (circles), 0.25 (up 
triangles), 0.088 (down triangles). 

In Fig. H we present two measured features of the lad- 
ders. In contrast to the case of a ID parallel array, there 
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are no additional lobes between / = and / = 1. Also, 
the critical current Iq remains relatively large, despite a 
/?l as low as 0.088. Similarly to Ref. [l^ we numerically 
solved the Eqs. (@)-(|) by using the same parameters as 
in the experiment (except for a, see above). In Fig. 
we compare the numerical simulations (solid lines) with 
the experimental data (symbols), which show good agree- 
ment. The /3l used in the simulations was calculated 
from the critical current of a single junction measured in 
the experiment. The calculations show some flattening 
at / = 0.5 for low /3l, which is also present in the experi- 
mental data. We found in simulations that in this region 
the ladder gets first filled with flux and only subsequently 
undergoes the depinning. We have observed this behav- 
ior only for low inductance, in good agreement with the 
analytical prediction of Ref. |ll| (that neglects inductance 
terms). In experiments with /3l = 3 and /3l = 0.88 we 
note the simultaneous presence of two different states at 
the same frustration value. We suppose that this is due to 
distinctly different initial conditions that can be realized 
in the ladder, while sweeping the bias current Ib through 
zero. This contradicts the prediction of Ref. [ll] where it 
is stated that the depinning of the whole ladder from a 
state different from the empty ground state occurs only 
around / = 0.5. The reason of this disagreement might 
be the neglect of the inductance in their calculations. 
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FIG. 3. Experimental (symbols) and numerical (solid lines) 
Jc vs / patterns of one of the ladders measured at different 
temperatures, in order to vary the critical current. The tem- 
perature has been derived from the gap voltage. Parameters 
are: N = 11; /3l = 3.0 (squares), 2.6 (circles), 2.1 (up trian- 
gles), 1.8 (down triangles). 

In the ladder with the largest /3l parameter (/3l = 3.0), 
we have measured the Jc vs / dependence also as a func- 
tion of the temperature T. At higher temperatures the 
decreased critical current causes a decrease of /3l- The 
results are shown in Fig. |3[ in physical units to underline 
the actual change of the critical current. Also in this case 



the agreement between the model and the experiments is 
rather good. 



IV. DISCUSSION 

We have characterized the static properties of Joseph- 
son ladders for different values of the self-inductance pa- 
rameter The experimentally observed dependencies 
of the critical current on frustration are in good agree- 
ment with the numerical simulations and show that the 
behavior of the ladder is clearly different from that of 
the ID parallel array (see Fig. |J). As it is well-known, in 
ID parallel arrays the critical (depinning) current is de- 
termined by the parameter /3l, and in the limit of small 
/3l the minimum frustration-dependent critical current 
is very small. Instead, the ladder critical current, even 
in the case of small /3l, never goes to zero. As it was 
already pointed out in Ref. [l^, with respect to ID paral- 
lel arrays, the presence of the horizontal junctions in the 
ladder leads to an "effective" increased /3£ ff , which for 
small discreteness can be by up to two orders of magni- 
tude larger than the natural fti, of the system, calculated 
(similar to ID arrays) from the junction critical current 
and the cell inductance. 

In order to show the particular mapping between 
Josephson ladders and ID parallel arrays, we carry out a 
simple quantitative analysis of the Eqs. (||). Let us con- 
sider the static case, when all Josephson junction phases 
are independent of time and satisfy the system of non- 
linear equations: 



?7(sin ipi 



■sinipi) = sinc/>i-7, i = 2,N. (4) 



By making use of the particular assumption that the hor- 
izontal phases tpi are small, we can eliminate the phases 
ipi from all equations and write the system of equations 
for phases <pi in the form: 



V 



[&_i - 20, + <p i+1 ] + 7 



(5) 



This system of equations coincides with the one describ- 
ing the static properties of ID parallel array (with hori- 
zontal junctions replaced by superconducting electrodes). 
The difference between a ladder and ID parallel array is 
that for the ladder we have now to use an effective pa- 
rameter 



0L + 2/77. 



(6) 



The deviation of /3£ ff from /3l originates from an addi- 
tional shielding (and vortex pinning) due to the presence 
of horizontal junctions, i.e. the horizontal junctions can 
accommodate part of the phase change. Thus, we ex- 
pect that this deviation disappears in anisotropic ladders 
when the critical current of horizontal junctions i^ or is 
much larger than the critical current of vertical junction 
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FIG. 4. Jc vs / dependencies of the ID parallel array with 
/3l = 2.7 (solid circles) and the ladder with /?l = 0.88 (open 
circles) . 

To verify the mapping given by Eq.(^|), in Fig. ^ 
we compare the patterns of a ID parallel array with 
/3 10 = 2.7 and an isotropic ladder with /3l = 0.88. For 
the ladder we expect /3£ ff « (3^ . The agreement is par- 
ticularly good at low frustration, but not in the vicinity of 
/ = 0.5, where the most critical assumption of our the- 
ory, i.e. is small values of the horizontal junction phases, 
breaks down. 

We would like to note here that a similar analysis for 
anisotropic ladders in the limit of small /3l has been car- 
ried out in Ref. 18, Moreover, for the case of ladder 



with three junctions per cell the mapping takes the form: 
/3£ ff = /?l + l/ 7 /- This mapping is in good accord with 
previously published data on the fa(/) dependence for 
ladder with three junctions per celo. We want to stress 
out here that this mapping is supposed to work only for 
the static case. In the dynamic state, when the Josephson 
vortices propagate in the ladder, the phases of the hor- 
izontal junctions start to oscillate and Eqs. arc not 
valid anymore, especially in the regime of large vortex 
velocities. Theoretical and experimental investigation of 
vortex propagation in Josephson ladders will be reported 
elsewhere£9. 



V. CONCLUSION 

We have reported measurements of the critical current 
in ladders of Josephson junctions. The /3l parameter has 
been varied by changing both the geometrical inductance 
and the critical current of the junctions. The results are 
in good agreement with numerical simulations, and show 
a behavior clearly distinct from the case of the ID parallel 
Josephson junction array without junctions in the hori- 
zontal branches. Using a simple quantitative analysis, we 
have shown that the static properties of ID parallel ar- 
rays and ladders can be mapped by properly scaling the 



sclf-inductance parameter This analysis well agrees 
with experimental data. 
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